Impurity induced spin texture in quantum critical 2D antiferromagnets 
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We describe the uniform and staggered magnetization distributions around a vacancy in a quan- 
tum critical two-dimensional S — 1/2 antiferromagnet. The distributions are delocalized across the 
entire sample with a universal functional form arising from an uncompensated Berry phase. The 
numerical results, obtained using quantum Monte Carlo simulations of the Heisenberg model on bi- 
layer lattices with up to » 10 5 spins, are in good agreement with the proposed scaling structure. We 
determine the exponent n' = 0.40 ± 0.02, which governs both the staggered and uniform magnetic 
structure away from the impurity and also controls the impurity spin dynamics. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Cx, 75.40.Mg 



Some of the most interesting physics of strongly inter- 
acting quantum systems arises in their response to im- 
purities. Metallic systems exhibit the Kondo effect, and 
much rich physics has been discovered in their response 
to a magnetic impurity which carries a localized spin 5. 
In contrast, Mott insulators have a particularly rich re- 
sponse to non-magnetic impurities. For example, the 
spin-gap compound CuGe03, which consists of dimerized 
pairs of 5 — 1/2 Cu ions locked into 5 = valence bonds, 
acquires magnetic order upon replacing a very small den- 
sity of the Cu with spinless Zn ions [l[ ; it is believed that 
an unpaired Cu spin is localized in the vicinity of the Zn 
impurity, and behaves like a localized 5 = 1/2 moment. 
The cuprate superconductors have also seen a variety of 
studies 0, H HH of the spin and charge correlations in 
the vicinity of Zn ions replacing the Cu ions within a su- 
perconducting layer; here there is also an unpaired spin, 
but its spatial distribution and dynamics are not fully 
understood. 

This paper will explore the effect of a non-magnetic 
impurity on a quantum critical Mott insulator, at the 
boundary between a magnetically ordered and a spin-gap 
state. We will describe the fate of the spatial magnetic 
structure of the localized impurity in the spin-gap state, 
as this state is tuned to the quantum critical point. We 
find that the strongly-interacting gapless excitations in 
the bulk lead to a non-trivial and universal spatial form, 
with power-law decay of spin correlations away from the 
impurity (analogous to skyrmions induced by dopants in 
the Neel state |6|). We will present new numerical and 
analytic results on the spatial spin distribution, building 
upon the scaling structure proposed in an earlier field- 
theoretical analysis @, [1| ■ Our results are of relevance to 
Mott insulators which can be tuned across the quantum 
phase transition. They also shed light on the cuprates, 
which are are in the vicinity of a magnetic ordering quan- 
tum phase transition. 

In our numerical investigations of a model spin-gap 
Mott insulator, we consider the spin- 1/2 Heisenberg an- 



tiferromagnet on a bilayer lattice. It is defined by the 
Hamiltonian 

H = J^S lli 'Si, j + J x ^S lii -S iti) (1) 

(*,i> » 

where S aj i is a spin-1/2 operator at site i on layer a = 1, 2, 
and denotes a pair of nearest-neighbor sites on an 
L x L open-boundary square lattice. With intralayer in- 
teractions in only one of the layers (a Kondo lattice) , as 
shown in Fig. [TJ the model has a quantum critical point 
when the ratio g = J±/J is g c = 1.3888(1) [9], with 
the spin-gap state present for g > g c , and an antifer- 
romagnetically ordered state for g < g c . There is con- 
vincing evidence 0, f]~oj | that this quantum critical point 
is described by the Wilson-Fisher fixed point of ip 4 field 
theory with 0(3) symmetry in 3 spacetime dimensions. 
Here, our discussion is conveniently presented in terms of 
the fixed-length formulation of this field theory, which is 
the 0(3) non- linear sigma model of the field unit length 
field n(r, r) (n 2 = 1), representing the local orientation 
of the antiferromagnetic order, with action 

5 bulk = ^J dr J d 2 r [(cVn) 2 + c 2 (V r n) 2 ] , (2) 

where r is imaginary time, r is the spatial co-ordinate, c 
is the spin-wave velocity, and the coupling g is a mono- 
tonic function of g. At the quantum critical point at 
g = g c , the correlations of the antiferromagnetic order 
are characterized by the power-law decay 

(n(r, t) • n(0, 0)} ~ (r 2 + cV)"^/ 2 , (3) 

with the exponent 77 rs 0.04 [ll[ of the Wilson-Fisher 
fixed point. 

Let us now add a non-magnetic impurity to the above 
systems. For the lattice Hamiltonian H, we create a va- 
cancy by removing a single spin. In a finite size Lx L sys- 
tem, this leads to a twofold degenerate 5=1/2 ground 
state. We take L odd and place the vacancy at the cen- 
ter of the top layer, as shown in Fig. Q] We choose the 
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FIG. 1: (Color online) An L — 5 incomplete bilayer model 
with a vacancy. The single unpaired central spin (yellow) 
constitutes frame R = 0. Frames R = 1 and 7? max = 2 
consist, respectively, of the surrounding red and blue sites. 

ground states to be eigenstates of the total S z ; in either 
of these states, all the (S z i t 2,i) are non-zero even in zero 
applied field (which we assume throughout), and exhibit 
an interesting spatial structure that we will describe. 

Next, we add the impurity to the field theory. The un- 
paired spin near the origin leads to a net uncompensated 
Berry phase between the antiferromagnetically oriented 
spins, adding an impurity term to the action [8[ 

Simp =isf drA[n(0, r)] • (4) 

which depends only upon the orientation of the order 
parameter at r = 0. Here S = 1/2 is the unpaired spin 
associated with the non-magnetic impurity, and A is the 
Dirac monopole function in spin space with V n x A = n. 
It was argued in Refs. @, Hf that Sbuik + Simp univer- 
sally describes the quantum impurity near g — g c . In the 
infinite system, the critical spin correlations are charac- 
terized by a different "boundary" exponent at r = 0: 

(n(0,r)-n(0,0)>~|r|-^. (5) 

A crucial property of the field theory <Sbuik + <S; mp 
is that, like the lattice model, for g > g c the ground 
state has total spin S. This means that for the con- 
served Noether magnetization density Q z (r), associated 
with the symmetry of 0(3) rotations of the action, has 
a non- vanishing expectation value even in zero field, and 
obeys 

(J d 2 rQ z (r)^ = S (6) 

in the ground state with maximum spin projection in the 
z direction. While this spin is localized in the spin-gap 
state with g > g c , there is a transition to a delocalized 
critical state at g — g c in which (as we describe below) 
the spatial extent of the magnetization is set only by 
the system size. Consequently, in an infinite system at 
g = g c we have {Q z (r)) = at all r even though Eq. ([6|) 
is obeyed [lit ]. 

A proper analysis of the delocalized spin texture in 
the ground state at g — g c requires imposition of a finite 



size L. We found from our numerical results, described 
below, that the spin distribution quite accurately obeys 
the scaling form 

(Q z (r)) = L- 2 $ Q (r/L), (7) 

where $q(jj) is a universal function (with no arbitrary 
scale factors) obeying J d 2 y$Q(y) — S. The same scaling 
form also emerges from a renormalization group analy- 
sis of the field theory 5buik + Simp, along with explicit 
results for $q(j/), and these will be described elsewhere. 
More useful here is the behavior of the scaling function as 
y — > 0, which describes the spin distribution in the vicin- 
ity of the impurity. A key feature of the theory [7] is that 
both the uniform and staggered spin operators of the bulk 
theory transmute into the same "boundary" operator as 
they approach the impurity. Here we are using n(0, t) 
to represent this boundary operator, and so the scaling 
dimensions implicit in Eqs. (0 [7]) suggest the operator 
product expansion lim r _^o Q( r , T ) ~ V \~ 2+v ^ 2 n(0, r), 
which in turn implies that 

<M^0)~zT 2 +"'/ 2 . (8) 

Similarly, we can also examine the distribution of stag- 
gered spin density, which is encoded in the spatial distri- 
bution of the order parameter n. From Eq. ([3|) we deduce 
the scaling form 

(n z (r)) =L-^/ 2 <$> n {r/L). (9) 

Now the operator product expansion to the same bound- 
ary operator is lim r ^o n ( r , T ) ~ |^| ^7 )/ 2 n(0, r), 
and this leads to 

$ n (y^0)~y-( 1+ "-^/ 2 . (10) 

Unlike <I>q, the integral of $ n is not quantized, and its 
overall scale is non-universal. Note that the theoreti- 
cal results in Eqs. (|7I8I9I10P are tightly constrained, de- 
pendent only upon a single exponent 77' in addition to 
the standard bulk critical correlation function exponent 
77 k 0.04. The value of 7/ has previously been estimated 
in the time domain, using Eq. (O [14| . We will show be- 
low that the numerics confirm that the spatial structure 
is also governed by this exponent, which we evaluate to 
higher precision than previously. 

In order to numerically study how the impurity- 
induced total magnetization S z = ^ S? = ±1/2 is dis- 
tributed in the system at T = 0, the lattice is decomposed 
into "frames" surrounding the vacancy, as illustrated in 
Fig. [TJ For each frame R we determine the uniform and 
staggered magnetizations, respectively, defined by, 



M (R) = 




MAR) = (^(-ir+^^-Sf A (12) 
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FIG. 2: The uniform frame magnetization Mq(R) (a) and the 
corresponding integrated quantity Io(R) (b) plotted according 
to the predicted scaling laws for several system sizes. In the 
inset of (a) we show the magnetization at R — 0, from which 
we obtain the exponent r/ = 0.40 ± 0.02. All other lines 
are shown using this value for the exponent. In (b) the data 
collapse is so tight that the L = 257 points almost completely 
hide the data for smaller L. 



where s = 2S Z = ±1 is included in order to make the 
contributions positive for both S z = ±1/2 states. Ex- 
pectation values are calculated in T > quantum Monte 
Carlo simulations utilizing the stochastic series expan- 
sion algorithm [l5| . Ground state results are obtained by 
choosing a sufficiently low temperature for each L. For 
the largest lattice we have studeied, L = 257, tempera- 
tures as low as T/J « 10 -3 are required for satisfactory 
T — > convergence ■ 

The uniform frame magnetization is related to the 
magnetization density Q z by Mq(R) oc RQ z (R) for 
R > 0. Hence, according to the scaling forms (|7|8|) . we 
should have for R/L small 
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(13) 
(14) 



where Io(R) is the integrated frame magnetization, 

R 



I (R) =^M (r), 



(15) 



which has to be exactly 1/2 at the edge of the lattice, 
i.e., for -Rmax = (L — l)/2. The uniform magnetization 
results are shown versus R/L in log-log plots in Fig. [5] 
The upper panel shows the frame magnetization, which 
is seen to collapse onto a single curve for L > 17, with the 
exception of the R = 1 points which scale with a differ- 
ent prefactor. Power-law behavior is seen for small R/L, 
including also the R = 1 data which approach such be- 
havior for larger L. For a fixed R, Iq(R) decays according 
to Eq. ([Til) as L' r >'l' 1 . This should hold also for R = 0, 
because /o(i? ma x) = 1/2. The inset of the upper panel 
of Fig. [2] shows the scaling of the R = magnetization, 
which in fact gives us the most precise determination of 
the exponent; rf = 0.40 ± 0.02. We use this value of 
7/ to draw the lines shown in the main panels of Fig. [5] 
The integrated magnetization Iq(R), shown in the lower 
panel, scales even better than Mq, with also data for the 
smallest lattices falling on the same curve. 

The theory is expected to capture accurately the be- 
havior for large lattices, far from the impurity and the 
edges, i.e., for large R but small R/L. It is therefore quite 
remarkable that even the Mq{R — 2) data fall on the 
common scaling curve and that the integrated magneti- 
zation scales almost perfectly even for very small lattices. 
The asymptotic power-law behavior is closely approached 
below R/L w 0.02. 

Next we consider the staggered component (fl"2"|) of the 
impurity induced texture and the corresponding integral 



I n {R)=J2 M n( r )- 



r=0 



Eqs. (PrO]) imply for small R/L 

Mk{R)l- {1 -^' 2 - 

I 7T (R)L-^/ 2 r. 



(3+v'-v)/2 



(16) 

(17) 
(18) 



r=0 



Results are shown in Figs. [3] We observe good data col- 
lapse of Mir for all L, and the slope vs R/L agrees very 
well with the value rj' = 0.40 determined above from 
the uniform magnetization. In this case the integrated 
quantity shows substantial subleading size corrections for 
small R, but an L — > 00 approach to the power-law be- 
havior shown by the line appears very plausible. 

In conclusion, we have presented analytical and nu- 
merical results for the uniform and staggered compo- 
nents of the spin texture induced by a static vacancy 
in a 2D quantum antiferromagnet at its quantum crit- 
ical point. The theory predicts scaling functions with 
asymptotic power-law behaviors, which are very well re- 
produced by the numerics. We have determined the value 
7/ = 0.40 ± 0.02 for the single exponent governing the 
asymptotic behavior of both the uniform and staggered 
structure. This exponent characterizes the influence of 



4 



1.0 



0.1 



- 

- 


1 






- 






■ 


7 — O 

/j — y 








□ 


L=17 


o 






o 


L=33 


8 - 









L=65 


o : 






• 


L=129 






(a) 




o 

111 


L=257 





0.01 



0.10 



R/L 



0.100 r 



0.010 



0.001 




0.01 



0.10 



R/L 



FIG. 3: Size and distance scaling of the staggered frame mag- 
netization Mtt (R) (a) and the corresponding integrated quan- 
tity Ik(R) (b). The over-all L dependence, L-f 1 "'))/ 2 an d 
£-(3->7)/2 £ or an( j respectively, has been divided out. 
The lines in (a) and (b) have slopes 0.68 and 1.68, respec- 
tively, corresponding to the value of the exponent jj' = 0.40 
extracted from the Mq (0) results in Fig. [2] 



the Berry phase in iSi mp and so, unlike the bulk the- 
ory, cannot be related to the exponent of any classical 
theory in one higher dimension. Our results support a 
central property of the "boundary" critical field theory 
0: the bulk operators for the staggered and uniform 
magnetizations transmute into the same boundary spin 
operator as they approach the impurity. Earlier time do- 
main studies [3] yielded r\' = 0.37 ± 0.05 [3]; the good 
agreement between the two approaches provides strong 
evidence that a single exponent indeed governs both the 
temporal and spatial impurity effects [7J. In addition to 
extracting the asymptotic power-law, our numerical cal- 
culations also give the full scaling functions for arbitrary 
distance from the impurity. We note that the integrated 
effects of the impurity are much stronger at criticality 
than in the symmetry-broken Neel state, where the in- 
duced disturbance of the magnetic structure around the 
impurity decays asymptotically as r~ 3 fl7j ]. 

NMR @, [3] and STM 0, i] experiments have probed 
the magnetization distributions around an impurity. Al- 
though these systems are not at a quantum critical point, 
we hope similar experiments in related system will be in 
the quantum critical regime, allowing measurements of 
the exponent rf using some of the observables discussed 
in Ref. 0. 
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